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Abstract 

In the context of softly broken Af = 2 super symmetric quantum chromo dynamics (SQCD), with a 
hierarchical gauge symmetry breaking SU(N + 1) — 1 -> U(N) — ^ 1, v\ 3> 1)2, we construct monopole- 
vortex complex soliton-like solutions and examine their properties. They represent the minimum of the 
static energy under the constraint that the monopole and antimonopole positions sitting at the extremes 
of the vortex are kept fixed. They interpolate the 't Hooft-Polyakov-like regular monopole solution near 
the monopole centers and a vortex solution far from them and in between. The main result, obtained in 
the theory with Nf = N equal- mass flavors, is concerned with the existence of exact orientational CP N ~ 1 
zero modes, arising from the exact color-flavor diagonal SU(N)c+f global symmetry. The "unbroken" 
subgroup SU(N) C SU(N + 1) with which the naive notion of non-Abelian monopoles and the related 
difficulties were associated, is explicitly broken at low energies. The monopole transforms nevertheless 
according to the fundamental representation of a new exact, unbroken SU(N) symmetry group, as does 
the vortex attached to it. We argue that this explains the origin of the dual non-Abelian gauge symmetry. 
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1 Introduction 



The last several years have witnessed a remarkable progress in our understanding of vortex configurations 
in spontaneously broken gauge theories, which carry continuous, non-Abelian internal zero modes: the 
non-Abelian vortices [H[2]. Their rich group-theoretical and dynamical properties have been subject of 
intense study [3]-[TT]. 

The physics of non-Abelian vortices is deeply related to the understanding of the concept of the non- 
Abelian monopole and that of the quark confinement, see e.g. Ref.[3j. Indeed, a detailed, fully quantum- 
mechanical analysis of Ad gauge theories with N = 2 supersymmetry has given important hints about 
the low-energy, effective dual gauge symmetry. In particular, fully quantum- mechanical light non-Abelian 
monopoles appear as the infrared degrees of freedom in the so-called r vacua of N = 2 supersymmetric 
QCD with Nf quark multiplets, playing the role of the order parameter for confinement (of non-Abelian 
variety) and for dynamical symmetry breaking |18j . 

The discovery of the non-Abelian vortex was partly motivated [2] by the desire to understand the 
physics of the quantum r- vacua, r = 2, . . . ,Nf/2, from a more familiar semi-classical viewpoint. In fact, 
semi-classically the connection between the vortex solutions and regular 't Hooft-Polyakov monopoles 
arises from the consideration of a hierarchical gauge symmetry breaking, e.g. 

SU(N + 1) U(N) 1 , vi » v 2 » A SU{N ) ■ (1-1) 

The monopole is supported by tt2(SU(N + 1)/U(N)) = Z; the low-energy vortex solutions correspond to 
non-trivial elements of iri(U(N)) = Z. The exact sequence of homotopy groups relates the two solitons 
of different co-dimension [4], and the global symmetry consideration tells us that the non-Abelian vortex 
implies non-Abelian monopoles sitting at its extremes. For oo > vi ^> v 2 > 0, one is inevitably led to 
consider metastable monopole- vortex complex solitons. 

The aim of this paper is to pursue further the study of such a monopole-vortex complex, including 
the numerical analysis of the field configurations involving both the magnetic monopole region and the 
asymptotic vortex-like region, with all fields approaching smoothly their vacuum expectation value (VEV) 
away from the complex. In this sense this paper is a continuation of the work by Auzzi et. al. [5]. We 
clarify also some aspects of the non-Abelian orientational moduli, extensively studied in the last several 
years in the context of the vortex solutions, and show how the properties of the non-Abelian orientational 
moduli can be extended to the whole monopole-vortex complex. 

The organization of the paper will be the following. In Sec. 2 after a brief introduction of the model 
we will focus our attention to the main characters of our paper: the vortex and the monopole, a particular 
care will be given to the study of the symmetries present at different energy scales. In Sec. 3 these two 
objects will be glued together and we will study the monopole-vortex complex as a whole. We will present 
both numerical and analytical results concerning the profile functions for the various fields and for the 
magnetic field, furthermore through a careful identification of the symmetries possessed by the complex 
we will be able to analyze the low energy dynamics of the orientational zero-modes living on the complex. 
The conclusions and results obtained in this work are presented in Sec. 4 while we give a more detailed 
analysis of the full set of equations of motion and the consistency of our Ansatz in the two appendices. 

2 The model 

To be concrete, we shall work with the softly broken N = 2 supersymmetric models, due to the many 
advantages they offer. The fields are those of the M = 2 gauge multiplet (the gauge field, the gauge 
fermion, the adjoint scalar and fermion) together with hypermultiplets. To be precise, we take an SU(N C ) 
(N c = N + 1) gauge theory with Nf = N flavors of hypermultiplets (quarks), and the mass parameters 
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are tuned such that at two hierarchically different scales the gauge symmetry is broken as in Eq. (jl.ip . 
The Lagrangian of the underlying SU(N + 1) theory has the structure 



£ = — Imr 

4"7T 



d^Tt ($t e -2V $ ) + 



d 2 9Tr(W Q W a 



+ £ (quarks) + f tfQ^Q* ^ ^ ^ 



£ (quarks) = ^ j £ Q (gtg-2V Q i + Q.JV + J #q ^Q.$Q< + m . QiQ 1 



(2.3) 



where mi are the bare masses of the quark fields and the complex coupling constant reads 

9 4vri 

The mass of the adjoint chiral multiplet \i breaks supersymmetry to M = 1. 

After elimination of the auxiliary fields the bosonic Lagrangian takes the form 
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[ mi + V2^q^ 
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2 

j 
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where j4 = 1,2,..., (iV + l) 2 — 1 and the sum over i = 1,2,..., Nf as well as a, b = 1, 2, . . . , N + 1 is 
implicit. In the construction of the monopole- vortex complex soliton solutions it turns out to be sufficient 
to consider the VEVs and fluctuations around them which satisfy 



[0 + ,0] = , 



i ~\i 



(2.7) 



therefore V\ can be set identically to zero in what follows. 

The vacuum expectation values (VEVs) of the scalar fields are determined from the minima of the 
potential following from Eq. (j2.6j) . e.g. see Ref. |18j . They are found to be 



% 
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for i = 1,2,..., r, a = 1, 2, . . . , N + 1 ; 



for i = r + 1 
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where the integer 



0,l,...,min{JV>,JV c -l}, N C = N + 1, 



(2.8) 
(2.9) 

(2.10) 
(2.11) 

(2.12) 



labels the possible (classical) vacua. The vacua of a given r are further classified according to which set 
of r (out of Nj) masses are used to construct the solution, leading to the total of 



#vacua 



3m{Nf,N c -l} 



r=0 



(2.13) 
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As explained in Ref. [18], by choosing a generic set of bare masses m; and by deforming with the 
M = 1 mass term ^Tr<3? 2 the continuous vacuum degeneracy is lifted altogether, leaving this discrete set 
of vacua. At small rrii and \x (<C A) the interactions become strong in the infrared, in all these vacua. 
These are indeed the vacua we are interested in0 

By tuning the bare squark masses mj to an equal, common value m, we see that an exact color-flavor 
diagonal SU(r) symmetry survives in a vacuum with a given r. For definiteness, below we shall work with 
the classical r = N vacuum where 

(0)= ncw -n ) ■ {q) = {q)=v2 {o lxN ) ■ ^ 

and 

Tfl 

v x = — = , v 2 = ^{N + l)mn . (2.15) 
v2 

are obtained by taking such a limit. Note that in this vacuum the gauge group SU(N + 1) is completely 
broken; at the same time, however, the color-flavor diagonal global SU (N) symmetry 

V U = ( U x ) oU^ , (^,40 = ( U j) x ) , (2.16) 

(U G SU(N) C SU(N + 1)) remains unbroken by both VEVs. It is an exact global symmetry of the 
whole system. The system is in the so-called color-flavor locked phase. A hierarchical symmetry breaking 
pattern (jl.ip is realized if 

|m|>H>A, .-. \vi\ > \v 2 \ ■ (2.17) 



Remarks 

(i) The terms containing the adjoint scalar mass \x play two crucial roles in our model. On the one hand, 

they induce the small squark condensates, (|2.14j) . bringing the system into a completely Higgsed 
phase. The existence of the vortex solutions in the low-energy approximation and their properties, 
all rely on this parameter. Note that due to supersymmetry, the high-energy approximate monopole 
solution (v2 = 0) and low-energy approximate vortex solutions (v% = oo) are both BPS-saturated. 

On the other hand, non-vanishing jj, introduces terms in V<i which make both the low-energy vortex 
and high-energy monopole "solutions" unstable (non-BPS). It is these terms which allow the two 
solitons of different codimensions to get smoothly combined into a monopole-vortex complex. 

(ii) Of course such a complex "soliton" is not a true solution of the field's equations of motion; it is only 

so under the condition that the monopole center positions are kept fixed. Under the assumption 
of a hierarchical gauge symmetry breaking (|2.17p . this is not a problem: it is a perfectly sensible 
(Born-Oppenheimer) procedure, as the motion of the massive monopole can be neglected in the first 
instance, in the study of low-energy fluctuations of orientational zero modes of the whole complex. 

(iii) Actually the case for working with non-BPS objects as these complex solitons can be made stronger. 
Just as in the case of the real- world quark-antiquark-chromoelectric string composites (the mesons) , 

1 This is one of the motivations for considering the system with generic masses and with /i / first, and then eventually 
taking the equal-mass or massless limit. On the contrary, if we considered directly the massless theory, or equal mass cases, 
we would find flat directions (continuum vacuum degeneracy); at a generic point along such Higgs branches, the coupling 
constant remain small at all scales. 
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our monopole-vortex-antimonopole complex will get stabilized after the quantization of the radial or 
rotational motions are appropriately taken into account. Of course, only the ground state for each 
flavor quantum number will be truly stable, as the pair production of the monopole-antimonopole 
from the vacuum, though suppressed, cannot be rigorously set to zero. But the same holds for the 
real-world meson states! The fact that these are not topologically stable as static configurations, is 
therefore not a shortcoming at all; rather it is the necessary price to pay to understand the real-world 
confinement mechanism. 

(iv) Our semi-classical consideration is valid for the parameter region 

|mj| = \m\ 3> S> A . 

As the bare quark masses and the adjoint scalar mass /i are decreased to small values of the order 
of A, the system becomes strongly coupled and the classical r = N vacua we started with turn into 
the weakly coupled, dual quantum r = vacua [18]. This means that the vortex (and monopole) 
orientational zero modes fluctuate strongly and Abelianize. A possible route to reach the quantum 
r vacua (2 < r < Nf/2), where light non-Abelian monopoles appear as the infrared degrees of 
freedom, through a careful tuning of bare masses and through the consequent vortex solutions with 
orientational moduli living on CP 7 * -1 x CP N ~ r ~ 1 , was discussed by some of us in Ref. [6]. 



2.1 Vortex solution far from the monopole center 



In the vacuum (|2.14j) the (N + l)-th color component of the squark fields has a mass of the order of vi 
and can be integrated away, in the study the low-energy dynamics. At mass scales much lower than v\ 
the theory is an SU(N) x U(l) gauge theory with Nf = N massless flavors, in the color-flavor locked 
phase. The vortex solutions in these systems have been a subject of an intense study for the last several 
years [3] -[17]. A lowest- winding vortex solution oriented in the (1, 1) direction in color- flavor spaed! has 
the form 



Qi(p) 



1 = ei i~jfi 



V\ + 




+ 



1 - /na(p) 



f N - 1 



y/2N{N + 1) 



N 



+ 



Ana(p) 



-1jv-i 
'N- 1 



1,2 



-N , 



^/2N(N - 1) 



LJV-l 



(2.18) 



with the profile functions satisfying the appropriate boundary conditions 

91,2(00) = V2 , /(oo) = /na(oo) = A (00) = Ana(oo) = , 

gi(0) = , d p 92(0) = d p A(0) = dp Ana(0) = , /(0) = / NA (0) 



1 . 



(2.19) 



We have introduced above the cylindrical coordinates 

P 



z 



\J x 1 + y 2 , tp = tan" 



-l (V_ 

x 



2 These vortices are precisely those studied by [7] for the special case of N = 2. The (jV + l)-th color component of the 
(massive) squark fields is vanishing and not shown in Eq. (|2.18[1 . 
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which are convenient in the description of the vortex configuration; the gauge field above in fact contains 
only the <p component, A v = — -[■■•] where [■■■] is the expression in the square bracket appearing in 
Eq. i£~E3- 

The n dependent terms in Eq. (|2.6p . coming from the adjoint scalar mass term /xTr<£ 2 , imply that the 
field <p behaves non-trivially (Eq. (|2.18p ): the vortex solution is necessarily non-BPS [7]. Its properties 
have been carefully studied by Auzzi et. al. [7] in the case of the low-energy SU(2) x U(l), N* = 2, theory. 
It has been shown that, independent of the sign of p,, or of m\ = rri2, the vortex tension is less than the 
BPS vortex (for p = 0). The behavior of the profile functions is quite similar, for small p, to the case of 
the BPS vortex analyzed earlier, except for the presence of the small non-vanishing profile functions for 
4> in and around the vortex core. 

For our purpose below, it is necessary to consider these vortex solutions in a gauge where the scalar 
quark fields do not wind at infinity. The SU(N) X U(l) C SU(N + 1) gauge transformation needed is 



-up 



j-y-(singular) 



IjV-l 



the vortex solution is indeed transformed into the form 

'n (p) 




/na(p) 
N 



f N-l 



-ljv-i 



0/ 



(2.20) 

(2.21) 
(2.22) 



while the <ft field remains invariant as in Eq. (|2.24j) . In this (singular) gauge all the topological features 
are hidden into the gauge-field singularity along the vortex core (p = 0), 




The magnetic flux is nevertheless regular everywhere and given by 

'1 



(2.23) 



1 d(pAp) 1 

Bz = a = ~ 

p dp p 



df(p) 1 
dp N 



ljV-l 



+ 



-N, 



dp N 



'N 



(2.24) 



with the B x ,B y components vanishing, so that the total flux is given by 



J~ : - 2ir J dp p B 




(2.25) 



Note that such a vortex leaves an SU(N — 1) x U(l) subgroup of the color-flavor diagonal global symmetry 
SU(N)c+f intact, and as a result develops orientational zero modes living in 

N-i SU(N) 



CP 



(2.26) 



SU{N-1) x 17(1) ' 

they are sort of Nambu-Goldstone modes, propagating along the vortex string. More concretely, there is 
a continuous infinity of degenerate vortex solutions related by 



JJ 



U 



1 



■u 



-1 



A 



u 



U 



1 



A; 



IJ- 



1 



U 



1 



u- 



1 



(2.27) 
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Note that the color-flavor N x N SU(N) matrix U belongs to the coset Eq. (|2.26p as an individual vortex 
solution (|2.22p preserves the subgroup SU(N — 1) x U(l) generated by 



fO \ / e i(N-l)a o 

VO e ia * TA ) ' V e~ ia °l 

where the T Al s are the standard SU(N — 1) generators. 

The existence of the exact orientational zero modes of Eq. (|2.27p possessed by the low-energy vortex 
solutions characterizes these solitons as non-Abelian vortices pQ 12] • They possess a genuine moduli 
space of solutions, all having the same tension. Their low-energy fluctuations (i.e. fluctuations carrying 
energies much lower than the characteristic vortex scale, v 2 ) can be shown to be effectively described as 
a two-dimensional sigma model on CP^ 1 = SU(N)/ SU(N — 1) x U(l). Dynamical features of these 
fluctuations, detailed moduli-space structures, and their group-theoretic properties have been the focus 
of considerable attention in the last several years [3]-[T7]. 



C U 



(2.28) 



Remarks 

(i) Note that even though the light squark fields have only the first N color components, the gauge field 

has a non- vanishing (N+l, N+l) element (see Eq. (|2T23j) or Eq. ([232]) ) due to the fact that the U(l) 
gauge group descends from the underlying SU(N + 1) gauge group. This turns out to be crucial 
in the consideration of the vortex-monopole complex below. In this respect, the present model 
(i.e. the model of [21 E] ) , where the SU(N) x U(l) theory arises as a low-energy approximation of 
a spontaneously broken SU(N + 1) theory, differs essentially from the genuine U(N) model studied 
by other groups PQ,[5]-[TT], where all scalar and gauge field components live in N x N color space. 
Even though they share many interesting features of the non-Abelian vortex solutions with our 
model, it is not possible in these latter models to relate the non-Abelian orientational moduli of the 
low-energy vortex to the notion of the non-Abelian monopole, which "lives" in the larger SU(N + 1) 
gauge group space. 

(ii) When the SU(N) and U(l) coupling constants are set to be equal - this would be the case for the 

theory just below the higher gauge-symmetry breaking mass scale v\ - the vortex solution above 
reduces exactly to the Abelian ANO vortex [?], embedded in the (1, 1) corner of the color-flavor space. 
We shall use below such a simplification e.g. for solving the monopole- vortex complex numerically 
(Subsection 13. 2D . but the more general discussion on the orientational moduli (Subsection I3.4p of 
the complex is independent of it. 



2.2 The monopole 

As the underlying SU (N+l) gauge group is simply connected, the vortex solutions reviewed above cannot 
be stable in the full theory, which contains the massive monopole excitations. A vortex must end at the 
two endpoints where monopoles related to the symmetry breaking SU(N + 1) — > SU(N) x U(l) are 
situated. The properties of such a monopole solution can be studied to first approximation by neglecting 
the small VEV, v 2 , or equivalently, by studying the system sufficiently close to the monopole center, i.e. at 
the distances R 

^->P~^, (2.29) 
\v 2 \ \Vl\ 

from the center. For q = = and for real <j) A neither V± nor V2 contribute, so the field configuration 
inside a sphere of such radius R should be approximately equal to the standard 't Hooft-Polyakov BPS 
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monopole solution [15] . embedded in an appropriate corner of SU(N + 1) gauge group. By choosing an 
577(2) C SU(N + 1) group generated by 



Si 















0/ 







Oiv-i 
/ 



'1 

Ojv-i 
,0 -1, 



(2.30) 



(which is broken to U(l) by the VEV of cf>) a monopole solution can be constructed explicitly as [20J: 

AT-l 



Mr) = A*{v) S a ; c/>(t) = (iV + 1) «i^ X (r) + v x 

r 



ljv-i 



2V-1 



(2.31) 



where 



Af(r) = e aji: ^a(r) , 



is the standard 't Hooft-Polyakov-BPS solution with 

gvir 



air) = 1 



smh(gvir) 
the latter behaving asymptotically as 

a(r) — > 1 , x( r ) ~^ sign(ui) 



x(r) = coth(^uir) 



gvir 



r — > oo 



(2.32) 



(2.33) 



(2.34) 



(for the antimonopole, x( r ) ~~ ^ ~x( r ))- The constant in the <j>(r) field is added so that it reduces asymp- 
totically to the vacuum expectation value of Eq. (|2.14|) . in a fixed (here chosen as (0,0, oo)) direction. 
The "magnetic flux" emanated from the magnetic monopole is given by 



n(s ■ r) 



(2.35) 



which of course is a well-known radial Dirac monopole field, embedded in the Si color directions. 

The "monopole" mass (the energy of the configuration around its center) can be approximately cal- 
culated a^l 



H 



Rewriting the Hamiltonian as 



d 6 xTt 



|r[<fl 



* (i^) 2 + 4^i 2 

2g z g z 



9 Z J\r\<R 



d 3 xTr 



-\Fij - e ijk (D k (f>)\ 2 + d k (e ijk Fij 



(2.36) 



(2.37) 



our monopole configuration is seen to satisfy approximately the non-Abelian Bogomol'nyi equations 

B k = V k <P , (2.38) 



3 The index a = 1,2,3 refers to the SU(2) group utilized to construct the solution; the gauge field At and the adjoint 
scalar (f> are both (N + 1) x (N + 1) matrices in the SU(N + 1) color group. 

4 We take the adjoint scalar field (f> to be real here and hence normalize it canonically as a real scalar field. 
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so that the monopole mass is given approximately by 

/ dS-Tr^B] = _(JV + l)| Vl , (2.39) 

J\r\=R 9 



*-4 

9 2 J|r|. 



where R is in the range (|2.29p . The monopole mass is, naturally, gauge independent. On the contrary, 
the "magnetic flux" (|2.35p is a gauge dependent quantity. To match the "vortex part" of the solution 
discussed in the previous subsection, it is necessary to choose a particular gauge in which the adjoint 
scalar field does not wind at infinity (the so-called singular gauge) in order to describe the complex. By 
introducing the spherical coordinates 

r 

- = (sin 6* cos sin 6* sin 99, cos6*) , 
r 

the required gauge transformation is given by 

<j> V <f> V* , At -»• V (Ai - idi) V ] , V = e -^ S3 e ieS2 e itpS3 , (2.40) 

where the Si's are the SU(2) C SU(N + 1) generators of Eq. (|2.30p . The adjoint scalar field is simply 
transformed to the fixed direction in the color space (chosen here in the 63 direction): 

jV-1 



{r) = (N + l)v 1 S 3X (r) + v 1 \ ljv-i , (2.41) 

AT-l 

whereas the gauge field has a slightly more complicated form. In view of our task below, of connecting 
the monopole solution to the low-energy vortex, it is convenient to express the gauge field in the new 
gauge in the components in the cylindrical coordinate system: 

A p = (£2 cos if — Si sin 99) (a(r) — 1) ; 

r 

A v = - [5 3 (1 -cos 9) -sin 0(Si cos ip + S 2 sin v?)(a(r) - 1)1 ; (2.42) 
P 

A z = —- (S2 cos ip — Si sin ip)(a(r) — 1) , 

r 

where the conversion between spherical and cylindrical coordinates involve the relation 

z = r cos 9 , p = r s'm9 , (2.43) 

while the azimuthal angle ip is common in the two coordinate systems. The gauge field displays the 
well-known Dirac string singularity along the negative z direction in this gauge: 

Wo=, - ~f ■ (2-44) 

The asymptotic (r S> t^t) behavior of the magnetic field B can be found by dropping terms multiplied 
by the factor (a(r) — 1) appearing in Eq. (|2.42|) : only A^ survives and 

B p ~ , ~ , B z ~ S 3 ^ , (2.45) 

which of course is a well-known radial Dirac monopole field, embedded in the S3 color direction. This 
radial magnetic field can be seen more easily by directly transforming the asymptotic magnetic field in 
Eq. (12351) by B { -> VB^, which yields 

* - ^ , (2-46) 
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The total magnetic flux through the surface of a sphere of radius R around the center, is then given by 



independently of R. 

3 Monopole-vortex complex 

The equality in the magnitude between the vortex flux (|2.25p through a plane perpendicular to the vortex 
axis, and the magnetic monopole flux (|2.47p through the surface of a tiny sphere around the monopole, is 
an example of the "flux matching" [U [12], showing that the vortex and monopole together form a smooth 
soliton complex. Below we elaborate this monopole-vortex complex in more detail. 

3.1 Generalities 

The first crucial observation is that such a smooth complex requires a compatible orientation between 
that of the vortex (in the color-flavor mixed space) and that of the monopoles (in the color space). For 
concreteness and for example the (1,1) corner of SU(N)c+f is selected for the winding for the vortex 
solution Eq. (|2.18p . and the associated monopole solution is accordingly embedded in the SU(2) subgroup 
of the color SU(N + 1) group, in the (1,N + 1) plane. Naturally, both the monopole and vortex can be 
rotated into any other directions by color-flavor transformations; in order to maintain the energy they 
must be rotated simultaneously (see below). 

The relative sign of the flux in Eqs. (12.250 and (12.470 indicates that the monopole studied in the 
preceding subsection is actually positioned at the rightmost endpoint of the vortex (as in Fig. |3J). The 
outgoing magnetic flux (in the direction of S s ) of the monopole is carried away in the vortex bundle on 
its left. 

Another crucial observation is about the behavior of the gauge field. The monopole field in the "sin- 
gular" gauge, Eq. (|2.42[) . exhibits the well-known Dirac string singularity along the negative z direction, 
Eq. (|2.44p . This behavior matches precisely that of the gauge field along the vortex core, Eq. f|2.23j> . 
confirming further the relative positioning of the vortex with respect to the monopole. The Dirac string 
singularity Eqs. (|2.23p . (j2.44j) is entirely harmless as it appears multiplied by the squark field in the squark 
kinetic term \T>^q\ 2 ; the latter drops to zero along the vortex core. The color magnetic field is smooth 
everywhere, as noted already. 

Note that the choice of the particular gauge (the so-called singular gauge) needed for matching the 
Dirac singularity of the gauge field, is precisely the one which guarantees the smooth field configuration ev- 
erywhere. In particular, in this gauge no fields "wind" around either the vortex axis or the monopoles. The 
scalar fields approach their vacuum expectation value everywhere outside the complex, as in Eqs. (|2.2ip 
and (J23H). See Fig. El 

In order to describe the field configurations interpolating the vortex and monopole, {q, A, <£} MV , we 




(2.47) 
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make an Ansatz of the form 

q 2 (p,z)l N -i 
cos 

A p = (52 cos if — Si simp) A(p, z) ; 

r 



N 



1 



In i + 

-N 



/na(p, z) 
N 



f N-l 



-ljV-i I — sin#(5i cos </3 + 52 sine/?) A(p, z) 




.4. 



sm( 



■(52 cos 9? — Si simp) A(p, z) ; 



<p= \Vl 



\(p,z) 



y/2N(N+l) 



(3.48) 



0/ 



The profile functions {q\,qi, f, In A, A, A, Ana}, with appropriate boundary conditions, can be determined 
numerically as we will do in the following section. 

3.2 Numerical solution 

In order to study these configurations numerically, we note that if the SU(N) and U(l) coupling constants 
are set to be equal our monopole-vortex complex is exactly the monopole-vortex complex generated by 
the symmetry breaking 



SU(2) -2* U(l) 



Vl > V 2 , 



(3.49) 



embedded in a larger color-flavor space. It is therefore sufficient for our purpose here to consider the 
minimal case, Eq. (]3,49p . The generators of the SU(2) group are t a = r a /2, where r a are the Pauli 
matrices and the scalar field has the form: 

where A a is the fluctuation around the VEV. The Lagrangian is then: 



1 



1 



where we have set 

£ = 4/um, v = 2V2~n, 
and our convention for the covariant derivatives and field strength is 



ml 2 - mr 6 + —^\ a T 
V2 



q' 
(3.50) 



T?a _ p, jta _ ft AO, _ abc Ab ac 

V„q = d^q + l -Alr a q ; 
V^ a = d^ a - e abc A b ^ c . 



(3.51) 
(3.52) 
(3.53) 



After the symmetry breaking at vi, the second color component of the squark field becomes massive, so 
we can set 

'qi(r,zy 
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The equations of motion for the system are: 



aba 



+ 



2 L 



q\T a V" 'q l - (V ' Qi y r 



£5 a3 + u\ a + q]r 



\ a i 



^= q jr a (ml 2 - mr 3 + ±=\K b ) J 



£5 a3 + uRe{X a ) + Tr (gV a g) ) (r a q) 



\>r a r 



, 1 

mlo - mr H =A a r 

V2 



(3.54) 
(3.55) 

(3.56) 



where we have defined \X\ 2 = X'X. In order to solve these equations numerically, we introduce an Ansatz 
adequate for the SU(2) theory, which is somewhat simpler than Eq. (|3,48p : 



A lr _ 



Z , ■ sf(P,z)-l 

2 ( r 2 cos f-nsmip) 



p 2 + z 

2 - 2 (n sin 93 - t 2 cos v?) 



p z + z z 



1 



(ri cos 99 + T2 sin 99) 



/(P^)-1 



1 

2p 



-V2m5 a3 + X a , X a = 5 a3 s(p,z) ; 

gi(r,z) 




(3.57) 



Substituting this Ansatz into the equations of motion, written extensively in Eqs. (jB.91|) - (|B.93p . one 
obtains the coupled differential equations Eqs. (|B.95ap -( B,95gp . A priori this system of equations seems 
to be an overdetermined system with respect to the function /(p, z). However as we have shown in the 
previous sections, the chosen Ansatz is well suited for both the monopole and the vortex, and hence we have 
assumed the existence of a solution to all the equations, solving only the system (|B.95dlB.95elB.95f|B.95g 



After the solution was obtained, we plugged it into the constraint equations (|B.96aHB.96bp as well as the 
remaining second order equation (IB.95b[) and indeed verified that the solution satisfies all equations. 

In solving the system of these differential equations the relaxation method is very useful. We introduce 
a fictitious time dependence into each of the profile functions and then write all the equations as: 



Ei 



dhj 
dt 



(3.58) 



where Ei denotes the equation of motion for the profile function hi, which is obeyed when Ei = 0. It is 
important that the equation Ei is of second order in spatial derivatives and that the sign of the Laplace 
operator is positive: Ei = d 2 hi + • • • with j summing over spatial dimensions. In this way the equation 
with the fictitious time-dependence resembles (a modified form of) the well-known heat equation, and 
when a stationary solution has been found, i.e. dt hi = the equation of motion for hi has been obtained. 

First we need to impose some reasonable initial conditions (i.e. a guess of the right solution) for the 
profile functions at t = and let the system evolve till a static solution is found. This was done by 
patching together three different solutions to the equations of motion. The first solution is the non-BPS 
solution of the vortex (of infinite length) , which is easily obtained numerically, see Fig. [2] and also [7] . 
The second solution needed is for the monopole and here it turned out to be most efficient to simply use 
the analytical solution of the BPS 't Hooft-Polyakov monopole. The third solution is simply the Higgs 
vacuum. Now patching the solutions together by using a piecewise function in the z, p coordinates, we 
obtained an initial "solution" , see Fig. [TJ 
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Higgs vacuum 
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C 




B 




B 






o 
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X 



z N non-BPS vortex solution N 

Figure 1: The figure shows how the initial conditions are made by patching approximate solutions together. 
A denotes the non-BPS vortex solution, B the BPS monopole solution and C is the Higgs vacuum. The 
boundary conditions are also shown, the thick line denotes the Higgs vacuum, the dashed denotes the 
vortex solution while the dotted line denotes Neumann boundary conditions. 



The final ingredient in the solving the vortex-monopole complex is probably the most crucial, namely 
the boundary conditions. Physical arguments tell us that the complex is not really semi-classically stable 
as it stands. Hence, one could worry that the relaxation method would shrink the system to merely the 
Higgs vacuum. In order to circumvent that, we implemented the boundary conditions for the system 
very similarly to what we did for Abelian soliton junctions in Ref. |21| . viz. fixing the monopole positions 
and hence the string length. In this way, the complex is a minimum of the energy - given the condition 
we impose. In Fig. Q]we have also shown the boundary conditions (BCs). The three sides far from the 
complex all have Dirichlet BCs, i.e. the Higgs vacuum. At the vortex string in the middle of the z-axis also 
Dirichlet BCs are imposed, namely the non-BPS vortex solution described above is used as the boundary 
condition - this is what fixes the string length. For the remaining pieces on the z-axis, we have imposed 
Neumann BCs due to cylindrical symmetry of the problem. 

The numerical calculation^] gave the results shown in Fig. [3] for the profile functions f(p,z), l(p,z), 
s(p,z), and q(p,z). Note that both scalar fields (the adjoint and fundamental) approach their VEVs 
smoothly far from the monopole-vortex complex in all directions. The monopole-vortex complex is im- 
mersed in the Higgs vacuum. Far from the monopole centers, the fields reduce to those of the (non-BPS) 
vortex, Fig. [2j whereas they reduce to the 't Hooft-Polyakov monopole solution in the small spherical 
region near the monopole center. 

These profile functions give rise to the magnetic field shown in Figs. [4] and [5j as expected the magnetic 
flux coming out of the monopole (isotropically near the monopole center and in the direction of S 3 ) is 
carried away in a vortex bundle located on its left, and gets reabsorbed by the antimonopole sitting at 
the far left extreme of the vortex. 

3.3 Macroscopic picture and duality 

Recently the monopole-vortex complex in the simplest gauge symmetry breaking ()3.49p 

SU{2) -^h U{\) Al, Vl > v 2 , (3.59) 
The numerical analysis has been done with the Mathematica package. 
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Figure 2: The profile functions in the transverse plane in the vortex region (far from the monopole). The 
solid line denotes the adjoint scalar field ^s(p), the dash-dotted line denotes the squark field q(p) and 
finally dashed line denotes the gauge field function £(p). The values of the parameters in this figure has 
been chosen as: g = l,m = 2,v = 2^2p = 0.1, £ = y/2~/5 ~ 0.28. 



has been discussed |22|, I23j in the London limit, in this limit the monopole is a point and the vortex a line, 
without width. In this approximation it is possible to perform the electromagnetic duality transformation 
explicitly, and the resulting dual theory can be solved for the electric and magnetic fields explicitly, in 
the presence of the monopole-vortex complex. In the presence of a static heavy monopole of unit charge 
sitting at r = 0, and with a semi-infinite vortex extending on its left, the magnetic (and electric if 9 ^ 0) 
fields are given by |23j 



E i = F °i = ^2 ^ (m ° n) ■ * = \^F jk = Bt° n) + B^5f , (3.60) 

where 

(mon) = l diG{r) Bivov) = rrf_ f° ^ ^ z _ z , 
9 9 J-oo 

and G(r) is the Green function, having the Yukawa form 

G(r) = — ^_ 2 5Hr) = e —, m= 9 -^ . (3.62) 
w -A + m 2 w r y/2 

Such a construction, in the case of a more general symmetry breaking 

SU(N + 1) SU(N) x U{1) 1 , 

(N > 2) and in the presence of a color-flavor diagonal symmetry, is not known; but as noted above, in 
the approximation g\ = the results, Eqs. (|3.60p - (|3.62p . represent the color-magnetic or electric flux in 
the 5 3 direction (Eq. (j2~3D"]) ). 

3.4 Orientational zero modes 

Let us now return to the symmetry breaking of our interest here 

SU(N + l) color ® SU(N) aavOT U(N) coloI ® SU{N) Ravoi SU(N) C+F , Vl > v 2 , (3.63) 
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(a) gauge profile function f(p, z) (b) gauge profile function £(p, z) 



2U 




-20 



(c) scalar profile function -^ s (p, z ) (d) squark profile function q(p,z) 

Figure 3: The behavior of the four profile functions for the gauge and scalar fields of Eq. (|3.57p . Note that 
the gauge and squark fields approach quickly the familiar vortex behavior: the gauge field is strongest 
along the vortex core where the squark field drops to zero. Note also that between the figures for £(p, z) 
and for q(p, z) the z axis is inverted in order to exhibit better their non-trivial (p, z) dependence. Away 
from the monopole- vortex complex all scalar fields reach quickly and uniformly their vacuum expectation 
values. The values of the parameters in this figure has been chosen as: g = l,m = 2,u = 2^2p = 0.1, £ = 
y/2/5 ~ 0.28. 
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Z 

(a) (b) 

Figure 4: The behavior of the magnetic field in the soliton complex near the monopole region. In (a) 
is shown a stream line plot while the intensity of the magnetic field is also shown in (b) by means of 
the color scheme. For negative values of the cylindrically radial coordinate p the plot is simply a mirror, 
i.e. in order to illustrate a cross section of the system. The values of the parameters in this figure has 
been chosen as: g = l,m = 2,v = 2^2/1 = 0.1, £ = ~ 0.28. 



V \ \ V ( - 



Figure 5: The magnetic field in the complete monopole-vortex-antimonopole soliton complex. For negative 
values of the cylindrically radial coordinate p the plot is simply a mirror, i.e. in order to illustrate a cross 
section of the system. The values of the parameters in this figure has been chosen as: g = l,m = 2, v = 
2y/2fi = 0.1, £ = y/2/5 ~ 0.28. 



15 



and consider the monopole- vortex complex (|3.48p . The color-flavor diagonal SU(N)c+f symmetry (I2.16P 
is an exact symmetry of our SU(N + 1) system in the Higgs phase: it is not an approximate symmetry 
of the low-energy SU(N) x U(l) gauge system. On the other hand, the monopole-vortex complex ()3.48p 
is not invariant under the full SU(N)q+f symmetry: it transforms as 



kU 

1 



U ^ ' /U " • (3.64) 



where U G SU(N)c+f- Of course, not all SU(N)c+f elements transform the configuration non-trivially. 
The complex (|3.48p is obviously invariant under the SU(N — 1) subgroup generated by 

A = 1,2, ...,(N-1) 2 -1 . (3.65) 




The U(l) subgroup of SU(N)c+f, generated by 

T( ° )= ( iY_1 -1 N ) ■ (3 - 66) 

acts in a subtler way on the monopole-vortex complex. The pure vortex configuration (|2.24p of the low- 
energy theory is clearly invariant under such a U(l) transformation, but the monopole-vortex complex is 
not obviously invariant. 

The problem is that the complex (13.48P lives in a larger, SU(N + 1) gauge space, i.e. it has necessarily 
non- vanishing color components in the (N + l)-th row and/or (N + l)-th column. However, decomposing 
Eq. (|3l)6j) as 



( T(0) )=T^ + (N-l)S 3 , T«=( 



N-l 



-1n-i I , (3.67) 

N-l 



in terms of two, (N+l) x (N+l) matrices, one sees that {q, A, (/>} MV is invariant under transformations 
whereas they transform under e lCiS3 simply as 

{q, A, </>} MV — > {q, A, </>} MV | {5i — > Si cos a — S2 sin a; S2 — > S2 cos a + Si sin a} . (3.68) 

This can be easily seen to be equivalent to the shift of the azimuthal angle 

ip — > if — a ; (3.69) 

in other words, it is equivalent to (i.e. can be undone by) a spatial rotation of angle a around the 
monopole-vortex axis, z. 

To summarize, non-trivial zero modes (|3.64p of the complex are generated by 

where the complex (iV — l)-component vector b represents the local inhomogeneous coordinates of 

N _i SU(N) 
CP ~ SU(N-l)xU(l) ' (3 ' 71) 

just as what was found by studying exclusively the low-energy vortex solutions (|2.26p - (|2.28p . 
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3.5 Demise and resurrection of the SU(N) symmetry 



In the small region around the monopole, (|2.29[) . the tiny VEV v% may be effectively set to zero; as the 
squark fields are then trivial (q = 0) the transformations ()3.64p look locally as global color transformations 
of the monopole solution. 

This shows that the isomorphism between the vortex moduli [2] and the monopole moduli! - both 
found to be CP^ -1 - is certainly not a coincidence. Nevertheless, attempts to understand the moduli 
space of the monopoles (representing iT2(SU(2)/U(l))) as something arising from its various possible 
embeddings in the larger color space G/H = SU(N + 1)/SU(N) x U(l), face inevitably the known 
difficulties: 

(i) Topological obstruction [23] : 

(ii) Non-normalizable gauge zero modes [25J; 

(iii) Non-local nature of Goddard-Nuyts-Olive (GNO) duality [26] . 

The first two issues have been discussed many times in the literature, and need not be reviewed here (see 
however below for observations concerning these points). It is perhaps worthwhile however to recall the 
third point, which is apparently the simplest of the three and at the same time, the deepest. The GNO 
quantization condition [26 

2/3-aGZ (3.72) 

where a are non- vanishing root vectors of H, and where the asymptotic gauge field is written as (in an 
appropriate gauge) 

Fij = UjkBk = e»jfc-|(/3 • H) , p = a* , (3.73) 

in terms of the Cartan subalgebra generators H of H, tells us that the set of degenerate monopoles 
are labeled by the dual weight vector f3. The duality implied by such a formula is clearly a natural 
generalization of the electromagnetic duality. The transformations of the monopole under the dual, 
magnetic group correspond to some non-local field transformations in the original description^. Simply, it 
is not the right question to ask how the monopole solutions associated with the gauge symmetry breaking 
G —> H transform under the "unbroken" global color subgroup H . 
Let us return to the symmetry breaking (|3.63j) 

SU(N + l) color ® SU(N) &avor ^> U{N) color ® SU(N) &aYOI ^ SU(N) C+F , ui > v 2 . (3.75) 

The "magnetic monopole" solution arising from the breaking at the scale v\ cannot be rotated without 
spending energy under the SU(N) co [ or : the latter is broken by the smaller squark VEV, t>2, Eq. (|2.14p . It 
costs energy to vary the embedding of the SU (2) subgroup Si within the SU (A r +l) co i or (Eq. (12.30f) ) . as the 

6 Here we neglect the U(l) zero modes related to the electric charge of the monopole; in the Higgs phase they survive only 
in the vicinity of the monopole center. See also Ref. [23] . 

7 The fact that in most papers in which these issues are discussed the gauge group is taken to be SU (N) appears to have 
helped to obscure this point, as the dual of SU(N)/Zn is again SU(N). It is however sufficient to generalize the discussion, 
for instance, to the case of a gauge symmetry breaking 

SO(2iV + 2) ( SO(2N) x 1/(1) , 

in this case the monopole transforms according to the dual of SO(2N) /1,2 ■ under Spin(2N). The monopoles are predicted 
to transform according to one of the 2 iv_1 dimensional spinor representations |13| . 
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squark fields have a preferred color-flavor locked direction (Eq. (|2.22p ). Forcing so (rotating the SU(2) co \ ov 
embedding) would distort the complex and yield a higher-energy configuration (an excited state). In other 
words, the subtle issues of the non-normalizable zero modes and the topological obstruction for defining 
the global color (i) and (ii), have been converted into an explicit breaking of the symmetry. 

However, the system has an exact color-flavor diagonal SU(N)c + p symmetry. Individual monopole- 
vortex complex configurations break it. They live in the continuous moduli of the solution space 

SU{N) C+F /U(N - 1) ~ CP"- 1 . 

The simultaneous color-flavor SU(N)c+f transformations acting on the original fields keep the energy 
unchanged; they induce movements on CP^ -1 . The complex transforms as in the fundamental multiplet 
N of an SU(N). 

This then represents a new, exact non-Abelian symmetry for the monopole: the latter transforms 
according to the fundamental representation of an SU(N) group. 

The fluctuations of the associated zero modes propagate only along the vortex; this makes sense 
because, according to the standard lore, the dual gauge group is in a confinement phase, the original local 
gauge group being in the Higgs phase. The monopoles are confined by the color-magnetic vortex, which 
in the dual picture represents the color-electric confining string. 

4 Conclusion 

In this paper we elaborated further on some earlier analyses [5j 0] of the monopole- vortex complex 
"soliton", arising in systems with hierarchically broken gauge symmetries. In a vacuum with an exact 
unbroken color-flavor symmetry (color- flavor locked phase), such a complex acquires orientational zero 
modes, analogous to those extensively studied recently in the context of non-Abelian vortices [3]-|17|. 
The field configurations are studied both qualitatively and numerically, verifying that they constitute a 
smooth "constrained" minimum of the energy, with the monopole and antimonopole positions kept fixed. 
The whole complex possesses the orientational degeneracy (arising from the exact color-flavor symmetry): 
it implies a new, exact continuous symmetry for the monopole, confined by the vortex. This, we argue, 
is the semi-classical origin of the dual gauge symmetry. Their fluctuations propagate in the monopole- 
vortex-antimonopole worldsheet strip, the monopole and antimonopole acting as the source or sink for 
the massless (Nambu-Goldstone like) excitations, propagating along the vortex in between. Extending 
the direct derivation of the 2c? vortex worldsheet sigma models on CP^ -1 [2j [9j \W\ (or in Hermitian 
symmetric spaces, SO(2N)/U(N), USp(2N)/U(N), etc., depending on the model considered |13j) to 
the case of monopole-vortex complex is a somewhat non-trivial task, presently under investigation. Our 
microscopic study here should serve as the basic starting point for such an analysis. 
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A BPS monopole and vortex solutions 

Here we will briefly check that the BPS equations for the monopole and the vortex are compatible with 
our second-order equations of motion. 
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A.l Monopole 

Let us start with the monopole; the BPS equation reads: 
while the covariant derivative of the latter can be written as 



i[(f>,T>j4>] , 



(A.76) 



(A.77) 



where we have used the BPS equation in the last equality. This equation agrees with Eq. (|3.54p in the 
BPS limit, with </> a real (taking into account a factor of 1/2 for the canonical normalization of the kinetic 
term) and setting q = 0. 

From the BPS equation (1A.76|) we obtain for the adjoint scalar 



2 ^ijk^k-Fij 







(A.78) 



due to the Bianchi identity. This equation agrees with Eq. (|3. 55 j) in the BPS limit; taking [i = and 
q = 0. 



A. 2 Vortex 

The vortex BPS equations are for SU(N) x U(l) 



fg =f e ij9 l [Tr(gVg) - £ 

Ffj =F e ijg 2 N TT{qh a q) = ; 
T>iq ± ieijVjq = , 



(A.79) 

(A.80) 
(A.81) 



where the upper (lower) sign gives a BPS vortex (anti-vortex) and a = 1, . . . , N 2 — 1 is the SU(N) adjoint 
index (in the N = 1 case only the first of the gauge equations appears). We consider here a symmetry 
breaking pattern like SU(N + 1) — > SU(N) X U(l) — > 1 and the BPS equations come about at low-energy 
in this system in the limit of fi — > while £ = — V2fi((f) } is kept constant. We will check this explicitly 
below. 



General case: SU(N + 1) -»• SU (N) x 17(1) -> 1 
The Lagrangian of the SU(N + 1) system is 



■4^(C) 2 + 7P^I 2 + IVI 2 -5 2 



1 



ml 



7V+1 



which leads to the equations of motion 

XyiW = /ABC L-Btp^C + ^B (p^t 



"«.Yf +A„i 



\/2^ A + «tt A ^ - V2g 2 qjt A (ml N+ i + v 7 ^) g 



V^q = -g 2 (V2fi Ke((f) A ) + Tr (qU A q)) (t 



1jv+i + V2(p 



(A.82) 

(A.83) 
(A.84) 
(A.85) 
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where A = 1, . . . , (N + l) 2 — 1 is the SU(N + 1) adjoint index and \X\ 2 = X'X. Let us now truncate 
the theory to the low-energy SU(N) x U(l) sector, i.e. we will integrate out the massive components of 
the squarks q and also the adjoint and gauge fields corresponding to broken generators, as they will have 
a mass of the order of \vi\. The vortex solution requires 4> a , a = 0, 1, ... , iV 2 — 1 (which is the adjoint 
field written in terms of the SU(N) x U(l) group with corresponding index a) to be fixed at the VEV 
(cfi a ) = 5 a0 (j) . As already mentioned the BPS limit corresponds to taking /i — > 0, with the product fj,(jp 
kept finite. We have: y/2(4> a )t a + m = 0. Hence, in this limit the second equation is trivially satisfied. 
For the first equation we have: 



VF a - 



-i 5 2 Tr [qU a V jq - (V jq )H a q] , 



(Ai 



We can combine the Eqs. (|A.79|) and (1A.80|) by setting the two coupling constants to be equal; go = gN = 9 
(and choosing the upper sign) which yields 



qq ] - 2^ C 







(A.87) 



Acting with the covariant derivative yields (note that it acts in the adjoint way on both terms) 



€ij 2 



(V iq )q^ + qV^ 



2 



{V jq )J - qV^ 



(A.88) 



where we have used the BPS equation (1A.81|1 in the last step. Multiplying by t a and taking the trace 
gives us Eq. (|A.86|) . 

The third equation reduces to: 



V(D iq = g 1 ( Tr(qH a q) - £d M ) (t a q) = 2- 



qq ] - 2^ C 



(A.89) 



with £ = -V2fi4>°. Using the BPS equation (|A~81~j) we obtain 



V(Diq = --eij[T>i,T>j]q = Fnq 



(A.90) 



Now upon insertion of the field strength (|A.87p into the above equation we see that it is exactly that of 
Eq. 



B Numerical analysis in SU(2) theory 

Writing out the equations of motion Eqs. (|3.54|) - (j3.56p explicitly in terms of the fundamental fields gives: 
Gauge fields 



d l diA 3a - did 3 A ia - e abc ( 2A ib diA jc + diA ib A jc - A b d j A ic ) + ( A ia A b A jb - A a ^A b A ib 



(B.91) 



„abc 



A ct d J A b - (d j X c ^)X b - V2m5 b3 (d j \ c + cPX 



2A ja \ b ^X b - A jb ( A a A 6t + A 6 A at 



+ y/2: 



m 



cS %ba ^jbS^ca xa3 xbc 



25 CA 8° 



5 a6 5° 



A ]b {\ c + A ct ) - 4m 2 {A> a - 5 a3 A^) 



r 

2 L 



q*&>q - (d j q^)q) + iA ja q*q 
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Scalar fields 



d' l di X a - 2e abc A ib di\ c - e abc d i A b \ c + A? A ib X b - \ a A ib A b - V2; 



m 



-t + q ] q) 5 a3 + v\ a ~\ - [b ab + ie abz ) X b q ; 



e ab3(Qi A bj + A ib A ib § a3 _ ^iS^ia 

(B.92) 



Squark fields 



d%q + % -{&J$)q + iA\d\ - \A™AU = (-£ + ^Re(A 3 ) + Jq) q - \ (|A 3 | 2 + \\ x + iA 2 | 2 ) 9 



+ + i(A\ + iA 2 )d\ = (Re(A 1 ) + iRe(A 2 )) g - (a^ + *A 2 t) g 



+ I(A3t( A i + iA 2 )-A 3 (A 1 t + a 2 t))g. 



(B.93) 



(B.94) 



Inserting the Ansatz (|3.57p into the equation rewritten in cylindrical coordinates we obtain the system of 
equations for the profile functions: 



+ P d P dj + vZ±z [(1 _ t) 8f + (/ _ 1)9i] _ 2 j^el + 4-4^/ = 

p p z + z z p z + z z 

z 2 (/-l)(l-^) 5 2 . o z(/-l) 3 z(/-l)(l-£) 2 w ^ n2 , 

- J}— > + y -z{f - l)q 2 + U + ^ ^ + 2z(f - 1 S - V2m 2 , B.95a 

P 2 ^/p 2 + z 2 2 p z + z 2 p 2 



pd;f + zd p dj + d p f - 



p 2 + z 2 



P 



: * [(/ - 1)9 Z * + (1 - I) d z f] - - ( ^-4 ) OJ 



P \p 2 + z 2 



z(f -l)(l-£) a 2 (t . 2i p(/-l) 3 , 1 



Pv 7 P 2 + z 2 - 



+ ^p(/-iy+ , , , 

' p z + z z p 



+ -{f - 1) (1 - i) 2 + 2p(/ - l)(V2m - s) 2 , (B.95b) 



p d p dj-zd 2 f+d z f--d p f+ 



z a ri 2z(f-l)d z £ , (/-l)(p 2 -z 2 )a/ 2(1-*)9 Z / , (p 2 + 2z 2 )(l-l)d p / 



77TI 2 



-+- 



pyJfF+Z 2 



■+- 



pa/p^T 

*(/-l)(*-2)£ 



, (B.95c) 



d 2 f + d 2 f + 2 



z(f-l)d p £ z(l-£)d p f ( P 2 -z 2 )d p f 2(f-l)d z £ , (!-£)£>,/ 2*0*/ 



+ 



PVP 2 + * 2 PvV + Z 2 P{P 2 + Z 2 ) ^/ p 2 + z 2 + z 2 p2 + z 2 



£(/ - itf + S (. - ^(/ - l) - ^-f V + ^ + ^ ' < B ' 9M > 

Z PVP +Z P 



Vp 2 + ^ 



p 2 + z 2 



2<lfl2( , Sp(f-i)(pd g f-zd p f) z(f-l) 2 (l-£)(f -l) 2 , 5 : 



9^ + 9 2 £ - -a^ + 
p p 



(p 2 + Z 2 ) 3 /2 



(p 2 + z 2 fl 2 



p 2 + z 2 



+ 



, (B.95e) 
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,2 ^ 2 1- 1 2 2 ^<?\ 2 t , , , (a - vgm)(/ - l) 2 

V + + -d p s = -g q s + — {q - f + i/s) + ^—^ 

p p 2 8 /r + z z 



* + + ^ = ^ + x^ 2 " * + + 2^T^ + • (B ' 95g) 



Note that plugging Eq. (|B.95cp into Eq. ()B.95ap yields Eq. (|B.95dp . 
Constraint equations 

2(/-l) 



pd 2 / - Z(9 P / H jJ- [pd z s - zd p s] = , (B.96a) 

s — \/2m 



1 



2(f - 1) [pd z q - zd p q] + q [pd z f - zd p f] + -(/ - + ^ " «)? = . (B.96b) 
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